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e-j . Abstract 



In the presence of primordial magnetic fields the oscillating cosmic axion field 
drives an oscillating electric field. The ensuing dissipation of axions is found to be 
inversely proportional to the conductivity of the primordial plasma. This coun- 
terintuitive result is essentially equivalent to "Zeno's paradox" or the "watched- 
pot effect" of quantum mechanics. It implies that the standard predictions of 
the cosmic axion density remain unaltered even if primordial magnetic fields are 
strong. 
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Besides neutralinos, axions are the only theoretically well-motivated particle can- 
didate for the ubiquitous cold dark matter that appears to be required in the standard 
picture of cosmic structure formation. Primordial axions were created by the "mis- 
alignment mechanism" |lj as well as by the relaxation of the string network which 
formed at the Peccei-Quinn phase transition ||. In units of the cosmic critical den- 
sity the relic axion abundance is found to be Q a h 2 = £ (1CT 5 eV/m) 1,175 , where h 
is the present-day Hubble parameter in units of 100 kms _1 Mpc _1 and m the axion 
mass. The exact value of the numerical coefficient £ = 0(1) is the subject of some 
debate 0. However, if axions are the dark matter in the galactic halos implied by 
astronomical observations, it appears safe to assume that their mass lies in the range 
10~ 5 eV < m < 10 -3 eV. The current round of direct search experiments || for the 
first time has a realistic chance of detecting galactic axions at the lower end of this 
mass range. 

Because of their nonthermal production, axions are essentially born as a Bose 
condensate, i.e. a classical, coherent field oscillation of the axion field. It is obviously 
important to understand if these oscillations are damped by dissipation effects which 
would thermalize, and thus reduce, the cosmic axion population. For example, it 
has been shown [|J that the thermalization by interactions with the cosmic plasma 
is inefficient for m < 10 _1 eV, corresponding to values of the Peccei-Quinn scale 
fa ^ 10 s GeV. (The axion mass and the Peccei-Quinn scale are related by m = 
0.62 eV 10 7 GeV / f a .) A proposed "coherent" damping mechanism || does not seem 
to be effective in practice |J. The possibility of resonant axion-photon conversion in a 
cosmological magnetic field of order 10~ 9 G has also been studied, with the conclusion 
that it yields no significant axion dissipation J7|. 

We presently study another dissipation mechanism which is expected if strong pri- 
mordial magnetic fields exist. They may arise during the early cosmic phase transitions 
[[|, and recently it has been shown that magnetic fields are indeed a stable feature 
of a second order (electroweak) phase transition ||. Locally the field could be very 
large. It is limited only by primordial nucleosynthesis arguments, which imply that 
B < 3 x 10 10 G at t ~ 10 4 s [fTU| . Because flux conservation implies that B ~ R~ 2 (R is 
the cosmic scale factor), at earlier times the field could have been much stronger. On 
dimensional grounds, a typical scale for magnetic field fluctuations should be B ~ T 2 
so that at the time of the electroweak phase transition local fields as high as 10 24 G 
could obtain. Depending on how such a large, random magnetic field scales at large 
distances, it could be the seed field needed to explain the observed galactic magnetic 
fields ||1 1|| . Let us remark that such large magnetic fields will not facilitate resonant 
axion-photon conversion in the early universe because the magnetic field removes the 
possibility for degeneracy in the refractive indices of the axion and photon fields. 

Magnetic fields would however couple to the cold axions and produce an oscillat- 



ing electric field, an effect which is used in the cavity experiments which search for 
galactic axion dark matter ||. In the early universe, such a field would give rise to 
a bulk velocity of the charge carriers (electrons, muons, and above the QCD phase 
transition temperature, quarks) and hence to a current. This induced current would 
rapidly dissipate by thermal collisions in the hot plasma. If the magnetic field is large 
this process might dissipate the cosmic axion energy density. Unexpectedly, however, 
we find this damping mechanism to be ineffective because the conductivity of the pri- 
mordial plasma is too large. We believe that the somewhat paradoxical nature of the 
axion survival story makes it worthwhile to communicate these results. 

In order to derive the equations of motion for axions coupled to the electromagnetic 
field with dissipation we start from the Lagrangian 

C = Id^a - \m 2 a 2 + ±E 2 - ±B 2 + g ai a E • B, (1) 
where a is the axion field, m its mass, 

_ a 

9a " = 2^Ta 

the axion-photon coupling constant, and f a the Peccei-Quinn scale. Because of the 
Nambu-Goldstone nature of axions, the Lagrangian Eq. ([!]) is valid only for a <C f a . 

The equations of motion for the coupled axion-photon system that follow from this 
Lagrangian have been derived by several authors fll2| . In our case they simplify if we 
assume that the electromagnetic field is dominated by a large homogeneous primordial 
magnetic field B and an induced electric field E while we neglect a higher-order induced 
magnetic field as well as the thermal radiation fields. If we include the possibility of 
an electric current density j the coupled equations of motion are found to be 

E = -g ai B a - j , 
a + m 2 a = g ai E ■ B . (3) 

Note that in the literature the terms proportional to g ai are often presented with an 
erroneous relative sign. 

The only conceivable macroscopic current density j is the one induced by the electric 
field E. Assuming a linear response of the medium we may use Ohm's law j = erE 
where a is the conductivity of the primordial plasma. From Eq. (§) it is then evident 
that B and E are parallel so that we may use E = |E| and B = |B| instead. This 
leads to our final equations of motion 

E = -Pa-aE, 
d + m 2 a = f3E, (4) 

where (3 = g a ^B. 



(2) 



The overall behavior of this system is perhaps easiest to understand if one uses the 
vector potential A as a dynamical electric field variable by virtue of E = —A. After a 
Fourier transformation the equation of motion is 




0, (5) 



\ — if3u —J 1 + m 2 
revealing that we have to do with two coupled harmonic oscillators of which one 



is damped, i.e. with an axion-photon mixing phenomenon |Tj|. This equation has 
solutions only if the determinant of the matrix vanishes, giving us the dispersion 
relations. 

One obvious solution is a static mode u)\ — 0, corresponding to a constant A and a 
vanishing a, i.e. to no electric or axion field at all. In the absence of dissipation (a = 0), 
there is a second static mode uj 2 = 0, even in the presence of mixing (/3 ^ 0). In the 
presence of dissipation this mode is purely damped, i.e. u>2 is always purely imaginary. 
For large dissipation (<r > m or (3) one finds uo 2 = —ia. The two remaining modes 
are = ±m if there is no magnetic field ((3 = 0), i.e. they correspond to the axion 
field. It obtains an electric field admixture for (3 7^ 0. 

In order to understand better the behavior of these mixed modes it is useful to 
consider a number of approximations. The natural oscillation frequency is the axion 
mass m. The "mixing energy" (3 is proportional to l/f a and thus to m; numerically 
j3 = 3.65 x 10~ 21 (B/G)m. We only consider magnetic field strengths small enough 
that always ^<mso that we are in the "weak mixing limit". Any deviation from 
0^4 = ±m will then be of order [3 2 . The shift of the real part of 0^4 is ±|/3 2 /m for 
cr = 0, and less for a > 0. The imaginary (damping) part of these frequencies is 

_ . , 1 f 2cr lm 2 for a < m, 

M«M = ¥ x , / f ' ( 6 ) 

[I/a for a ^> m. 

Therefore, if a -C m the axion field is damped more strongly for an increasing conduc- 
tivity as naively expected. In the "strong damping limit" a ^> m, on the other hand, 
the actual damping rate of the axion modes decreases with increasing a. 

Actually the conductivity of the primordial plasma is huge. In the regime m e -C 
T <^ Tq CD one finds for an isotropic relativistic electron gas fl4j| 



^plas T . . 

a = ~ , (7) 

where CG> plas is the plasma frequency and cr co n the collision cross section. This result 
is valid for fields smaller than the critical field B c = m 2 /e = 4.41 x 10 13 G, above 
which the electrons cannot be treated as free, and the conductivity Eq. (|7]) should be 
multiplied by a factor B / B c . 



In the nonrelativistic regime, relevant for the recombination time, one should use 
the conductivity of a nonrelativistic isotropic hydrogen plasma, given by 

(2T) 3/2 



" " K 3/2 " 1/2 - ( § ) 



Either way, the conductivity is very large compared with the axion mass so that we 
are always in the "strong damping limit." 

Because the axion field amplitude is dissipated away at the rate T = /3 2 /(2a), 
the rate for the dissipation of axion number density is that of the squared amplitude, 
T = ft 2 /a. This enables us to check whether this damping mechanism is effective or 
not in an expanding universe in the following way. 

The damping of axions would be effective only if the dissipation rate for the axion 
number density exceeds the Hubble expansion rate, i.e. if damping would occur on a 
time scale fast compared with the cosmic expansion time scale. Assuming magnetic 
flux conservation (3 scales as 1/R 2 and hence essentially as T 2 , while a scales as T in 
the relativistic epoch so that T oc T 3 . The expansion rate is roughly T 2 jnip\ (Planck 
mass mpi) so that axions must be dissipated early if they are dissipated at all. One 
crudely estimates T/H ~ (/3 /m) 2 m 2 mpi/T 3 . Even if f3/m were near unity at the QCD 
epoch, the factor m 2 mpi/Tq C - D is sufficiently below unity for the axion masses under 
consideration. Therefore, axions are not significantly damped by this mechanism. 

This conclusion remains valid if one takes the expansion of the universe explicitly 
into account, and also if one considers the nonrelativistic epoch as we shall now demon- 
strate. In order to study the explicit solutions of Eq. fl3|) it is useful to take 9 = a/ f a as 
a dynamical variable describing the axion field. With the initial conditions O(to) — 1 
and E(to) = 0, the approximate solution for /? C m C a is 

6(t) ~ exp(- /32(t ~ to) ) cos[m(t - t )}. (9) 

1(7 

However, this still neglects the expansion of the universe which we shall now include. 

Assuming a flat Robert son- Walker metric with R(t) the dimensionless scale factor 
of the universe, we define 

Foi = —REi, Fij = eijkBkR 2 . (10) 

Including dissipation as before, one may then derive the following equations of motion: 

E = - 9 -^-2HE-aE, 
R z 

e + 3Fe + m 2 e = ^§. (11) 

Here g = a/(2n) = g aj f a and H = R/R is the Hubble parameter, and we have made 
use of flux conservation which entails R 2 B = const = Bo, so that in a comoving volume 
the ratio of the magnetic energy density to radiation density remains fixed. 



In this case one cannot perform a simple mode analysis as in the nonexpanding case. 



Instead, eliminating E from Eq. (11), one obtains a third order differential equation 
for 0: 

+(7H + a)Q + (m 2 + (3 2 R~ 4 + 12H 2 + 3H + 3Ha)& + (AHm 2 + am 2 )& = (12) 

where we have used (3 = g aj B = gB / f a . Assuming a radiation dominated universe, 
we may find an approximate solution to Eq. ( |T2"D by writing 0(i) = C(t)Q (t), where 
0o (t) is the solution in the absence of expansion, and C(t) is a slowly varying function. 
Note that a ^> m ^> H except exactly at the QCD phase transition, above which the 
axion mass is actually smaller than H. The mass starts to turn on very rapidly, how- 
ever, so that the hierarchy of scales is valid when T < Tq CD . Writing a = a (t /tY^ 2 , 
to lowest order in {3 and H = 1/ (2t) we then obtain 



// \ 3/4 

9(i) ~ (J>) exp 



ffio ( /^0\l/2 
(To Vt ) 



cos[m(t — t )] • (13) 



One readily observes that damping cannot compete with the expansion of the uni- 
verse, which redshifts the energy of the background magnetic field, and at t — > oo 
damping tends to an asymptotic value. Thus the decrease in the axion field density 
due to magnetic field catalysed dissipation in the early universe is completely negligi- 
ble. For example, for a critical field B ~ 10 13 G one finds, substituting Eq. (|7|) into 
Eq. fll3|) , that dissipation reduces by a factor (1 — 5) x 10~ 22 . This result will remain 
qualitatively unchanged even if we allow for magnetic fields greater than B c . 

In the nonrelativistic regime the conductivity is given by Eq. (§) and H = 2/(3t), 
whence Eq. (|l3l) is replaced by 



0(0 ^(j) exp 



3/? 2 ^o / /to 



4a Vt 



!2)2/3 _ i 



cos[m(t - t )] • (14) 



Therefore, damping effects are again negligible. 

The standard predictions of the cosmic axion density remain unaltered even if there 
exist strong primordial magnetic fields. The reason is that the electric field never has 
a chance to grow large because the ohmic dissipation is so effective. The pump turning 
the axion energy into electric energy produces a mere trickle. 

This rather counterintuitive result is closely related to "Zeno's paradox" or the 
"watched-pot effect" JOJ. If two quantum states mix, the transition rate between 
them is suppressed if one of them is repeatedly measured, causing the system to 
remain "frozen" in the original state. In our case axions would oscillate into photons 
because the two fields are coupled by the external magnetic field. The dissipation of 
the electric field energy can be viewed on the quantum level as a photon absorption 
and thus as a "measurement" of the system to be in the electromagnetic state. If it is 
measured too frequently, it stays frozen in the initial axion state: the transition rate 
is suppressed inversely proportional to the rate of measurement. 
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